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Abstract 

Statistical theory of the dielectric susceptibility of polar liquid crys- 
tals is proposed. The molecules are not uniaxial but similar to cones. 
It is assumed that the permanent dipole moment of a molecule is par- 
allel to the axis of the rotational symmetry. The ordering of the phase 
is described by means of the mean field theory based on the Maier- 
Saupe approach. The theory was used to calculate the temperature 
dependence of the order parameters and the susceptibilities. Predic- 
tions of the model for different sets of parameters are investigated. 

PACS number(s): 61.30.Cz, 77.84.Nh 

1 Introduction 

Liquid crystal phases consist of anisotropic molecules, either elongated or 
disk-like [lj. In many theoretical [2j [3] and computer simulation studies 
[U [5] it is assumed that the molecules have the (or D 2 h) symmetry 

and that the interactions do not change if a molecule will point to the oppo- 
site direction. However, real molecules are always less symmetrical. Somoza 
and Tarazona [6] showed that the molecular asymmetry can influence the 
packing fraction, the relative density jump and the order parameter at the 
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isotropic-nematic transition for hard body systems. Recently Barmes et al 
[7] reported results from Monte Carlo simulations investigating mesophase 
formation in systems of hard pear-shaped particles. For different elongations 
isotropic, nematic, interdigitated smectic A 2 phases were obtained and some- 
times glassy behaviour was seen. A pear-shaped molecules was also used by 
Davidson and Mottram [8J in a continuum theory model of switching in a 
bistable menatic liquid crystal device. 

In the molecular-statistical theory of Maier and Saupe the potential en- 
ergy of interactions is 

V(ni,n 2 ) = v + v 2 P 2 {n x • n 2 ), (1) 

where Pj are the Legendre polynomials, the vectors ni,n 2 determine the 
orientation of molecules. There is a nematic phase for v 2 < (vq is not 
important). We would like to investigate the interactions that depend on the 
sense of the vectors ni,n 2 (it corresponds the symmetry of molecules) 

V(f£i, n 2 ) = v + «iPi(f?i ■ n 2 ) + w 2 P 2 (ra! ■ n 2 ). (2) 

The parameters Vi can be choosed by means of the excluded volume method 
described by Straley [9]. As a result we get always v\ > 0. The energy of 
interactions will be lower for n\ = —n 2 and this leads to the conclusion that 
only isotropic and nematic phases can be obtained. 

We can obtain negative (or positive) v\ and new ordering of the phase if 
we consider interactions of two electric dipols. Under the assumption that 
p = pn, the interactions have the form 

TT ^ p 2 [{ni-n 2 )-Kni-&){n 2 -K)} 
U 12 {u,m,n 2 ) = , (3) 

47T€ U 6 

where u = uA is the distance between dipols. We would like to get average 
interactions that do not depend on the vector A 

U^(u,nx,n 2 ) = — J dKUw&n^n^g&ni,^), (4) 

where the function g(A, f?i, n 2 ) describes the neighbourhood of the molecule 
1 (how often we can find the molecule 2 in the direction A at the constant 
u). The normalization is 

— JdKg(K,ni,n2) = l. (5) 
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For isotropic distribution (spherical molecules) g(A,fti,n 2 ) = 1 and we get 



U^(n 1 ,n 2 ) = 0. (6) 

But for anisotropic molecules it is more reasonable to try different form of 
the function 

g(A,n 1 ,n 2 ) = 1 - aP 2 (n 1 • A) - aP 2 (n 2 ■ A), (7) 

where a is related to the anisotropy of molecules: rodlike (disclike) molecules 
have a > (a < 0). 

U 12 (n u n 2 )- ^— -5 . (8) 

Thus, for disclike molecules (a < 0) one can use V\ < and this will 
produce ferroelectric ordering for sufficiently low temperatures. Generally 
we can say that v i describes the result of sum of the excluded volume effects 
and dipolar interactions. Our aim is to investigate the influence of the v± 
term on the ordering and the dielectric susceptibility of the phase. 

As far as the susceptibility is concerned, we would like to recall some 
results in order to provide the context for our considerations. The potential 
energy of the dipole p in the electric field E is U = —p ■ E. For a system of 
iV noninteracting dipole moments in the volume V we get the polarization 
P from the Boltzmann distribution 

P = ^L ((3 P E) , (9) 

where /3 = l/(k B T) and L(x) = cothrr — l/x is the Langevin function. For 
small values of the argument, L(x) ~ x/3, and thus 

P = (10) 

3V y 1 

The susceptibility is defined as 

dP Np 2 

eoX -dE-Wk^T- (U) 

We get a 1/T dependence called the Curie law. It indicates that the dipoles 
are noninteracting. 
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For the anisotropic polar molecules and anisotropic phases the suscepti- 
bility is a tensor 

eoX*p = Wfi . (12) 
The polarization can be written as 

P a = ^r(n a ), (13) 

where p{n a ) is the average value of the dipole component in the a-direction 
in the presence of the electric field. Hence two different averages are involved. 
The bar refers to the potential energy of the dipole moment in the electric 
field and the brakets to the nematic potential. We calculate it using the 
linear approximation 

(n a ) = (M 1 + PpnpEp)) = (3p(n a np)Ep. (14) 

The susceptibility is 

eoXa/3 = — y—{n a n p ). (15) 

The components of the susceptibility in the nematic phase oriented along the 
z axis are 

Np 2 (l + 2S) 

eo Xzz = wkaT , (16) 
A^ 2 (1 ~ S) 

€oXxx ~ 3Vk B T ' 1 ' 

where S = (^2(^2)) is the order parameter. Note that the average suscepti- 
bility is 

X = i(x„ + 2 XBW ) = ^ = x fa0 (18) 

as for the isotropic phase. The equality \ = X ts ° is correct for nonpolar 
molecules but for polar ones we know from the experiments that \ < X %so ■ 
More advanced calculations were given by Maier and Meier [TU] who extended 
the Onsager theory of the susceptibility to nematic liquid crystals. If we 
neglect the induced polarization we can write their results as 

(3x±l\NAl + 2S)_ 
60X22 \2x + 3j 3Vk B T ' [W) 

e QX * x - (20) 
0X \2x + 3 3Vk B T ' V ; 
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where x is the average susceptibility of the phase. 

The organization of this paper is as follows: In Sec. [2] the mean field 
theory of the nematic ordering is provided. In Sec. [3] the formulae for the 
susceptibility are derived. Section @] is devoted to some applications of the 
presented theory. Section [5] contains a summary. Appendixes provide the 
definitions and main properties of the basic functions (App. A) and the ex- 
cluded volumes for molecules of different shape (App. B). 

2 The mean field theory 

The potential energy of molecular interactions of the molecules is given by 

V(R 1 ,R 2 ) = £ VjESiR- 1 ^), (21) 
i=o,i,2 

where Ri are the sets of the three Euler angles describing the orientation of 
the molecule i, E^J are the basic functions defined in App. A. The molecular 
orientation in the phase is described by the distribution function 

JdRf(R) = l. (22) 

The mean of some function A = A(R) we calculate as 

(A) = J dRf(R)A(R). (23) 

The state of the system is described by order parameters (E$) but the most 
important is S — (Eq ). The internal energy of the system is 

N r 

U = - J dR 1 dR 2 f(R 1 )f(R 2 )V(R u R 2 ), (24) 

whereas the entropy of the system is written as 

S = -k B N J dRf(R) ln[f(R)C N }. (25) 

The energy of permanent dipole moments p = pn in the internal electric field 
E int is 

U E = N(p- E int ). (26) 
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Generally, the internal electric field differs from the external field E ext [TT] . 
Here we assume that both field are equal E mt = E ext = E. 
The total free energy of the system is the sum 



F tot = U + U E -TS. 
In the mean-field approximation we get a potential energy 

3 » 

and from the Boltzmann distribution we get 

f(R) = exp[-pW{R)]/Z 1 
where Z is a normalization constant. The consistency condition 



W{R X ) 
leads to equations 



dR 2 f(R 2 )V(R 1 ,R 2 )-p(R 1 )-E 



w 



U) 



(EjS)vj - p6 lj (EJ lfl - EyS-xu + E z 8 .) 



(27) 



(28) 



(29) 



(30) 



(31) 



It is useful to introduce the dimensionless parameters Sjj$ = — fiw^Q for j > 



ln/( J R) = EE^M. ) W- 

3 M 



(32) 



Sqq is responsible for the normalization and it depends on other S$ with 
j>0 



o(0) 
D 00 



In 



J>0 ii 



J Q0 



W 3* = ^ 

v 9S <$ 



{E^ E u0 ) - {E fJj0 ){E 1/0 ). 



The equations have the form 



(33) 
(34) 
(35) 

(36) 



The solution is orientationally stable only if the matrix 



l>0 P 



(37) 



is positive definite. The isotropic phase is orientationally stable if fiv\ > —3 
and (3v2 > —5. 



3 The dielectric susceptibility 

The dielectric susceptibility tensor is defined by the eq. f|T2|) and the orien- 
tational polarization we calculate as 

P a = (n a )Np/V. (38) 

Note that the polarization depends on the electric field via the distribution 
function. The components of the susceptibility are 



Np • ^dS^ 

j>0 M 



<W„ = -f (39) 



where the derivatives are calculated from the equations 

EE[M^ + «J]# = PpSvSor, (41) 



k>0 v 



EEMp + Wl^f = /W*- ( 4 2) 

fc>0 v UCj x 

Now we are in the position to disscuss the results for different possible phases. 
3.1 The isotropic phase 

For the zero field all order parameters are equal to zero. For the nonzero 
field the phase has the symmetry (the symmetry of the electric field). 
Thus we have to use the parameters Sqq only. 

Np 2 

eT = e oXzz = e oXxx = + ^ . (43) 
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For vi < we get the Curie- Weiss law describing the divergence of x when 
we approach the Curie temperature from above. For v± > the susceptibility 
is finite. 



3.2 The uniaxial nematic phase 

For the zero field the phase has the symmetry and the parameters Sqq 
with j even are present. For the parallel field the phase has the symmetry 
Coot, and the parameters Sqq with all j are present. For the perpendicular 
field the phase has the symmetry C 2v and we have the parameters Sqq with 
j even and with j, fi odd. The expressions Wj£ with j + k and /i + v 
even are nonzero. 

eoX = Nv2ml1 (44) 

eoXzz v(i + Vl w^y [ } 

oXxx ~ v(i + (3 Vl wuy [ } 

Note that for v i > we get x < X %s ° ■ 

3.3 The ferroelectric phase 

For the zero field and for the parallel field the phase has the symmetry 
and the parameters Sqq with all j are present. The expressions Wj£ with 
/i + v even are nonzero. 

oXzz { v ) (i + fawtfxi + pvtwft) - pvtfvtiwjgr {0> 

For convenience, we express the important elements Wf£ by means of the 
order parameters (Pj) in the case of the ferroelectric phase 





= [l + 2(P 2 )]/3-(P 1 ) 2 , 


(47) 


w 12 


= [2(P 1 ) + 3(P 3 )]/5-(P 1 )(P 2 ), 


(48) 


w% 


= [7+10(P 2 ) + 18(P 4 )]/35-(P 2 ) 2 , 


(49) 


wll 


= [1-W1/3, 


(50) 


w\l 


= [(Pi) - (P 3 )]v / 3/5, 


(51) 


wll 


= [7 + 5(P 2 )-12(P 4 )]/35, 


(52) 
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= o, 

= [7-10(P 2 ) + 3(P 4 >]/35, 
= 0. 



(53) 
(54) 
(55) 



4 Exemplary calculations 

In this section we carry out calculations for different physical systems of 
polar molecules described by the considered model. The phase diagram of 
the model is shown in Fig. 1. Three phases are present: isotropic, uniaxial 
nematic and ferroelectric. Dashed half-lines in the picture denote different 
physical systems with the fixed parameters V\ and v% (in fact, the fixed ratio 
t>2 to Vi is important). The parameters of the model were calculated by means 
of the excluded volume method. The formulae for molecules with different 
shapes are given in App. B. 

Let us start the description of the results from the system of uniaxial 
molecules similar to ellipsoids or cylinders [half-line (d) in Fig. 1, v± — 0]. 
In all pictures, T denotes the dimensionless temperature. T = 1 corresponds 
to the transition from the isotropic to the nematic or ferroelectric phase. 
The susceptibilities are scaled in order to obtain x = 1 for T = 1. The 
temperature dependence of the order parameter (P 2 ) and the susceptibilities 
are presented in Figs. 2 and 3, respectively. There is the first order transition 
from the isotropic to the nematic phase. On decreasing the temperature 
the susceptibility splits into Xxx < Xzz- Xzz runs to the infinity whereas 
Xxx remains finite. This is typical for the uniaxial molecules with the dipol 
moment parallel to the symmetry axis. 

Let us consider less symmetrical molecules similar to cones with H/D = 4, 
where H and D denote the height and the diameter of the base, respectively 
[half-line (e) in Fig. 1, t> 2 = —6%]. Note that only the isotropic and nematic 
phases are present. This reflects the fact that liquid crystals are in general 
nonpolar. On the other hand, one should expect the flexoelectric effect which 
is polar [13J. The temperature dependence of the order parameter (P 2 ) is sim- 
ilar to shown in Fig. 2. The temperature dependence of the susceptibilities 
is presented in Fig. 4. We have known splitting of the susceptibility but here 
both Xxx and Xzz are finite. In Fig. 5 the results for t> 2 = —v\ are given 
[half-line (f) in Fig. 1] where Xxx and Xzz are almost the same order. Note 
that this value of the t> 2 /t>i cannot be obtained from the excluded volume 
method but we use it in order to understand the behaviour of the model. We 
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conclude that if cones become longer then their properties are more similar to 
ellipsoids and cylinders. For shorter cones it is more difficult for the electric 
field to rotate the molecule and that is why the susceptibility \ zz is lower 
(but still greater then Xxx)- 

A ferroelectric phase is present in the considered model when the param- 
eter V\ is negative and there is a resultant tendency to arrange the dipols 
in a parallel fashion. The temperature dependence of the order parameters 
and the susceptibilities in the case of v?, = 5v\ [half-line (c) in Fig. 1] are 
presented in Figs. 6 and 7, respectively. There is the first order transition 
from the isotropic to the nematic phase at T = 1. Next, there is the second 
order transition to the ferroelectric phase at T = 0.763. 

The temperature dependence of the order parameters and the suscepti- 
bilities in the case of v-i = v\ [half-line (b) in Fig. 1] are presented in Figs. 8 
and 9, respectively. There is the first order transition from the isotropic to 
the ferroelectric phase. Xzz is always finite and it has the strong maximum 
at the transition. 

For the case of v 2 = we get the simple equation for Sqq only [half-line 
(a) in Fig. 1] 

S$ + V1 L(SS) = PpE z . (56) 

In the zero field there is the second order transition from the isotropic to the 
ferroelectric phase at Tc = — i>i/(3&b). The order parameters shown in Fig. 
10 can be expressed as 

(Pi) = L{S$>), (57) 
(P 2 > = 1-T/Tc. (58) 

As far as the susceptibility is concerned, we recover a general mean field 
result: the susceptibility in the ordered phase is half of that in the disordered 
phase (at the neighbourhood of the transition point) |14j . The reason is that 
the effect of an external field on the polarization must be smaller in the 
ordered phase, due to the existing internal field in that phase, which has a 
stabilizing effect on the polarization. In Fig. 11 we plot the inverse of the 
susceptibility as a function of the temperature. The dashed line shows the 
results from the Landau description for the ferroelectric phase. We have an 
expected agreement in the neighbourhood of the transition point. 
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5 Summary 



In this paper we considered the mean field theory with interactions which are 
not uniaxial. We showed that they correspond to molecules similar to cones 
or molecules with permanent dipol moments parallel to the axis of rotational 
symmetry. Three different phases were obtained: isotropic (I), nematic (N) 
and ferroelectric (F). We found the first (I-N, I-F) and the second order 
transitions (N-F, I-F) for some sets o parameters. 

We also derived the formulae for the dielectric susceptibility in the case 
of polar molecules. We neglected the induced polarization and the difference 
between the external electric field and the internal electric field acting on a 
molecule. The susceptibilities depend on the density, square of the dipolar 
moment, temperature and the order parameters. The anisotropy Ax = Xzz~ 
Xxx is positive as a result of the assumption on the dipol moment orientation. 

Our calculations show that for positive or small negative values of v\ the 
ordering is the same as in the case of v\ = (uniaxial molecules). Below 
the isotropic phase there is the nematic phase with the -Doo/i symmetry. On 
the other hand, the nonzero v\ strongly changes the component of the sus- 
ceptibility which is parallel to the symmetry axis of the molecule. We also 
find that for v\ > the average susceptibility is lower than the isotropic 
susceptibility what is in agreement with the experiments for polar molecules. 

The presented model describes only liquid crystal phases with transla- 
tional symmetry. For real liquid crystals below the nematic phase (for lower 
temperatures) one can usually find the smectic phase translationally ordered 
in one direction (rod-like molecules) or the columnar phase translationally 
ordered in two directions (disc-like molecules) [TJ]. Sometimes the nematic 
phase is even absent. It is possible to extend our model for smectic phases 
by extension of the McMillan model [16] or to include the columnar phases 
by extension of the model by Feldkamp et al [17]. Work is in progress on 
the theory were many limitations of the presented approach is removed, i.e. 
any orientation of the dipol moment is possible, the induced polarization and 
difference between external and internal electric fields are taken into account. 

Appendix A 

Below we list the properties of the functions Ej^J. The functions can be used 
to describe any physical quantity which depends on the three Euler angles. 
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1. The definition is 

-y \2+<5 /J+<5o 



Ej$(R) = [j 5 ) U " U l[(l + .) + (l-.)sign( /U )sign(^)] 
x[D$(R) + sign( / i)sign(z/)(-l)^ J D°') V (R) 



+ sLga(v)(-l)"D%}_ v (R) + sign( / i)(-l)^ (j '^( J R)], (59) 

where R = ((f), 9, if>) (the three Euler angles), j is a non-negative integer, 

fx and v are integers. Functions D^J are standard rotation matrix 
elements jTJ] and 

• f \ \ 1 for x > . . 

Slgn(x) = ( -1 for x < 0. (60) 

Note that 

sign(-x) = -sign(:r) + 25 0x . (61) 

2. The functions EjfJ are real. 

3. The functions satisfy the orthogonality relations 

iiv wp* yR) = S jk 5, p 5 ua 8n 2 /(2j + 1). (62) 

4. Let us assume that the three Euler angles R = (0, 9, i/j) determine the 
orientation of the three unit orthogonal vectors (l,m,n). The coordi- 
nates of the vectors can be expressed by means of the functions E^J 

l x = e[^ (R) = cos 9 cos cos ip — sin sin ip, 

l y = —-Bin (R) = cos $ s i n cos "0 + cos s i n 

Z z = E'oi^i?) = — sin 9 costp, 
m x = -E^! (R) = — cos 9 cos sin ip — sin cos ip, 
m y = E^l_i(R) = — cos 9 sin sin ip + cos cos -0, 
m z = i?Q_i(i?) = sin # sin -0, 
n x = Eiq(R) = sin # cos 0, 
n y = —E^^R) = sin 6* sin 0, 

n z = = cos#. (63) 
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5. The functions are proper for the description of uniaxial molecules. 

We can express them as products n a x r^ y rf z ^ where b = or b = 1. 



-^oo — - 1 ? 

-t^lO — n x, 

-^00 _ n z? 

F (1) - -71 

£/_io — %, 

4? = l/2(-l + n 2 2 + 2n :c 2 )v / 3, 

#10 = n x n z VS, 

E$ = -l/2 + 3/2n 2 2 , 

-E-io = —n z n y y/3, 

-£'-20 = n x n 2/V^3, 

= l/4^(-3 + 3n/ + 4n :c 2 )v / 5v / 2, 

4? = l/2n g (-l + n g 2 + 2n x 2 )V3V5, 

E$ = l/4n x (-l + 5n z 2 )V3V2, 

E$ = -3/2n z + 5/2n 2 3 , 

^So = -1/4^-1 + 5™/)^^, 

-^- 3 2o = n x n z n y V3V5, 

E% = -l/4{-l + n z 2 + 4n x 2 )n y V5V2, 

E { £ = 1/8(1 - 2n 2 2 -8n x 2 + n z A + 8n x W + 8n x 4 )v / 35, 

4? = l/4n x n z (-3 + 3n 2 + 4n x 2 )V35V2, 

E$ = 1/4(1 -8n 2 2 -2n x 2 + 7n 2 4 + 14n x 2 n 2 2 )v / 5, 

= l/4n :c n 2 (-3 + 7n 2 2 )v / 5v / 2, 

= 3/8-^V + fn/, 

= -l^ra^-S + Tra^v^A 

Si 4 i = l/2n x n y (-l + 7n z 2 )y/5, 

E% = -l/An z n y {-l+n 2 + 4n x 2 )v / 35v / 2, 

£i 4 Jo = l/2n x n y (-l + n z 2 + 2n ;E 2 )v / 35. 
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Appendix B 

Excluded volumes V exc i for uniaxial molecules with the D^h and sym- 
metry. V mo i denotes the volume of a molecule. The excluded volumes were 
used to calculate the parameters of the interactions Vj (j = 0, 1, 2) from the 
expressions 

V exc i{e z , e z ) is proportional to v o + f i + v 2, (64) 
V exc i(e z , —e z ) is proportional to v — v-y + v 2 , (65) 
V exc i(e z ,e x ) is proportional to vo — V2/2. (66) 

1. Ellipsoids with the axes 2a x 2a x 2c (Dooh symmetry). 

Kno* = (4/3)vro 2 c, (67) 
V ea;c /(e;,e;) = SV mo i, (68) 
We.,4) « (4/3)7r(a + c) 2 (2a). (69) 

2. Cylinders with the length H and the diameter 2R symmetry). 

Vmoi = kR 2 H, (70) 
V exd (e z ,e z ) = 8V mol , (71) 
V exd (e z ,e x ) = AttR 3 + 27tR 2 H + ARH 2 . (72) 

3. Cones with the height H and the diameter of the base 2R (C^ sym- 
metry) . 

V mol = (l/3)7rR 2 H, (73) 
V exc i(e z ,e z ) = UVmoh (74) 
V exc i{e z , -e z ) = 8V mo i, (75) 
V exd (e z ,e x ) « 12R 3 + 8R 2 H + ARH 2 . (76) 

4. " Towers" that consist of two cylinders which have the common axis of 
rotational symmetry (C^ symmetry). The cylinders are of the height 
H/2 and the diameters R and 2i?. 

Vmoi = (5/8)7rR 2 H, (77) 

Vcxcife, e z ) = 10V moh (78) 

V'exd^, -e z ) = Wmol, (79) 

V^d^, e x ) « 15R 3 + 12R 2 H + ARH 2 . (80) 
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Figure 1: Phase diagram of the model considered in the paper. Three phases 
are present: isotropic, uniaxial nematic and ferroelectric. Dashed half-lines 
denote different physical systems: (a) v 2 = 0, (b) v 2 = v±, (c) t> 2 = 5v\, (d) 
v\ = 0, (e) i>2 = — 6t>i and (f) v 2 = —v\. For a given physical system on 
decreasing the temperature we are moving from the center (0,0) to the edge 
of the figure. 
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Figure 2: Temperature dependence of the order parameter (P 2 ) for v± — 
[half-line (d) in Fig. 1]. There is first order transition from the isotropic to 
the uniaxial nematic phase. T denotes the dimensionless temperature. 
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Figure 3: Temperature dependence of the susceptibilities for v-y — [half-line 
(d) in Fig. 1]. On decreasing the temperature the susceptibility splits into 
Xxx < Xzz- For T — >■ Xzz becomes infinite. T denotes the dimensionless 
temperature. 
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Figure 4: Temperature dependence of the susceptibilities for t> 2 = — 6v± 
[half-line (e) in Fig. 1]. Both susceptibilities are finite. T denotes the 
dimensionless temperature. 



19 




Figure 5: Temperature dependence of the susceptibilities for t> 2 = —V\ [half- 
line (f) in Fig. 1]. This value of the i^/fi cannot be obtained from the 
excluded volume method. T denotes the dimensionless temperature. 
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Figure 6: Temperature dependence of the order parameter (Pi) and (P2) 
for t>2 = 5i>i [half-line (c) in Fig. 1]. There is the first order transition 
from the isotropic to the uniaxial nematic phase at T = 1 and next, the 
second order transition to the ferroelectric phase at T = 0.763. T denotes 
the dimensionless temperature. 
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Figure 7: Temperature dependence of the susceptibilities for vi = 5v\ [half- 
line (c) in Fig. 1]. In the nematic phase we have Xxx < Xzz- At the transition 
from the nematic to the ferroelectric phase Xzz diverges (the Curie- Weiss 
law). The vertical dashed line denotes the transition point between nematic 
and ferroelectric phases. T denotes the dimensionless temperature. 
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Figure 8: Temperature dependence of the order parameters (Pi) and (P 2 ) 
for t> 2 — v i [half-line (b) in Fig. 1]. There is the first order transition from the 
isotropic to the ferroelectric phase. T denotes the dimensionless temperature. 
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Figure 9: Temperature dependence of the susceptibilities for v 2 = v 1 [half- 
line (b) in Fig. 1]. Xzz is finite and it has the strong maximum at the 
transition from the isotropic to the ferroelectric phase. T denotes the dimen- 
sionless temperature. 
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Figure 10: Temperature dependence of the order parameters (Pi) and (P 2 ) 
for v 2 = [half-line (a) in Fig. 1]. There is the second order transition 
from the isotropic to the ferroelectric phase. T denotes the dimensionless 
temperature. 
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Figure 11: Temperature dependence of the inverse of the susceptibilities for 
V2 = [half-line (a) in Fig. 1] . The susceptibility Xzz in the ferroelectric phase 
is half of that in the isotropic phase (in the neighbourhood of the transition 
point). The dashed line shows the results from the Landau description. T 
denotes the dimensionless temperature. 
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